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_JL ■ The temperature dependence of the surface tension j(T) at liquid gal- 

lium is studied theoretically and experimentally using light scattering from 
capillary waves. The theoretical model based on the Gibbs thermodynamics 
relates dj/dTto the surface excess entropy density —AS. Although capil- 
lary waves contribute to the surface entropy with a positive sign, the effect 
of dipole layer on AS" is negative. Experimental data collected at a free Ga 

0^ ' surface in the temperature range from 30° C to 160° C show that d'y/dT 

changes sign near 100 °C. 
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The temperature dependence of the surface tension of liquids is very im- 
portant in applications such as Marangoni convection and crystal growth. 
However, little is known about the specific surface forces influencing this 
dependence in liquid metals. 

Our previous works [|], |2| were concerned with the surface tension of 
mercury as a function of temperature, which was studied using quasielastic 
light scattering from capillary waves ||]. Although the experimental data 
show a decrease in surface tension with T, the theoretical explanation of this 
fact seems to be non trivial. The recent discovery of surface layering in liquid 
gallium and mercury indicates an ordered metal structure perpendicular 
to the surface which is about a few atomic diameters thick [Q, (5J. This 
ordering can drastically influence the entropy density profile in the surface 
zone, reducing the entropy density at the surface compared to the entropy 
density of the bulk. It was suggested |I]] that a liquid metal surface is a 
two-component system comprising quasi- free surface electrons and positive 
ions. In this case d^/dT would depend on the surface excess entropy as well 
as the chemical potential of surface electrons M. However, the ionic and 
electronic profiles of any complicated shape should be approximated by the 
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same right-angled profile P] in order to preserve the electrical neutrality of 
metal as a whole. Thus, the two-component model is inapplicable to a free 
metal surface, whose excess entropy should be evaluated in the framework 
of a one-component model [||, f| 
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where AS is defined AS = ^ OQ dz{S{z) - S b 8(z)). Here 0(z)=O for z> 
and 0(z) = l for z< 0; z=0 denotes the surface position and the integration 
is performed from the liquid bulk (z=-oo) to a vapor phase (z=oo) and 
S(, denotes the entropy density in the liquid bulk. We split the surface 
excess entropy into three parts describing surface layering, capillary waves, 
and a surface double layer J2]. As expected the contribution from surface 
layering into AS is negative |1 1 , whilst that contribution from capillary waves 
is positive. It was suggested || that the electric field in the double layer 
would contribute to the surface excess entropy with a positive sign if the layer 
thickness varies with T but the surface charge density does not. The present 
paper demonstrates that this is not the case; the thickness of the dipole layer 
remains insensitive to the temperature variations and the contribution from 
electric field into AS" should be negative indicating a decrease in the surface 
electron density on heating. 

2 Theoretical estimates of AS 

A liquid metal comprises two components: free electrons and positive ions. 
Free electrons behave as a quantum medium: even at zero temperature their 
energy is nonzero, and is usually written via the Fermi energy Ej: E = -E/N, 
where N is the number of electrons. The spatial electronic density profile 
p e (z) has a nonzero width (even at OK), that leads to some redistribution 
of electrons between bulk and the surface. In other words, some electrons 
are ejected from the bulk and concentrate on the vapor side of the interface. 
This charge separation leads to a surface double layer, suggested by Frenkel 
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Unfortunately, no theory describing the surface tension of liquid metals 
is well established. Numerical simulations of the electronic and ionic den- 



sity profiles [11] predict oscillations of the density profile but the driving 
force of these oscillations is not very clear. Only one paper, to the author's 
knowledge, treats this problem analytically |l2jl . Due to its importance for 
the present context we repeat the main results of this paper. 

A detailed analysis done by Samojlovich at K ensures that the elec- 
trostatic force —p e V(j) leads to Maxwellian elastic stresses II 22 = —H xx = 
€qE 2 /2 (z axis is directed along the surface normal and x-along the sur- 
face) or electrostatic pressure P=€qE 2 /2 deforming the ionic fluid. This 



isotropic pressure acting on the topmost layer makes the atoms there over- 
compressed compared to atoms in the bulk. These considerations predict a 
non-monotonic density profile with a maximum in the topmost layer, which 
was demonstrated recently on the surface of liquid gallium and mercury. 

The surface tension of any liquid can be expressed via the tensor of 
anisotropic stresses with the components U zz and Yl xx appearing as a result 
of the density gradient in the surface zone. The surface tension of a liquid 
metal at OK is given by 

/co /-co 1 F) n 

dz(H zz - U xx ) = / dzlxv-i./) 2 - e E 2 } (2) 
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where % w = %— [12J. This expression is obtained from the balance between 
the surface forces, namely the electrostatic force —pS7(f> and the quantum 
pressure of electron gas given by von Weizsacker's term Xw J dz^(Vp) 2 . If 
only the electrostatic energy were accounted for the surface tension would be 
negative [jlj]. In order to calculate the surface tension one should know the 
electron density profile and the solution of Poisson's equation in the surface 
zone. 

The temperature behavior of the surface tension of a liquid metal is a 
most intriguing problem. While the changes in tension for nearly all liquid 
metals, as expected from literature, are only few percent over 100 K |lj|, the 
temperature derivative d'j/dT should be extremely sensitive to the density 
profile Q • In the following text we will show some basic relations connecting 
the surface excess entropy to other physical properties of the surface layer. 

It is supposed that the surface entropy comprises different parts J7]]: the 
first, pertinent to all fluids, is due to capillary waves. The second part 
is an orientational entropy AS or determined by surface layering. As we 
discussed in Q the variation of the entropy should have a minimum in the 
surface zone because the surface atoms are more ordered than atoms in the 
bulk. This would require AS or < as it was schematically drawn in [0]. 
There might be an orientational contribution in AS" due to ordering in the 
surface plane. This ordering, which would appear in the form of a hexatic 
monolayer on isotropic liquid substrate, is responsible for a change in sign 
of the temperature derivative of the surface tension of normal alkanes and 
some liquid crystals |jsj, [?J . For liquid metals in-plane surface order was never 
reported in experiments; only atomic ordering along the surface normal was 
detected |4], [5| . Also from a theory this hexatic layer on a pure liquid metal 
cannot be expected: for the formation of the Wigner lattice one needs a 
very low electron density |l(J. We show below that the impact of surface 
layering on AS is naturally linked to the specific surface effects namely the 
anisotropic pressure and the dipole layer. 

According to Frenkel's analysis the surface excess entropy is || 

AS = n s k B \ii^. (3) 



where n s is the number of atoms per unit area, O is the mean frequency 
of capillary waves, defined as ui c = — ^ max dqoj(q)2irq and ui p is the mean 
frequency of bulk phonons. A S is governed by the ratio of two mean fre- 
quencies surface and bulk and is positive because u) p > Q c . In order to 
calculate Q c one should know the dispersion relation oo(q) in the whole q- 
range: from long waves to the shortest wave limit. It should be borne in 
mind that the well-known relation uj ~ q 3 ' 2 derived from a linear theory of 
capillary waves is valid only for long waves, ie. relatively small q. For high 
wavenumbers, giving the main contribution in u> s , this relation is not valid 
any more because linear theory fails at some q such that the wavelength of 
the wave is of the same order as the amplitude. Hence it is not possible to 
calculate correctly the contribution of capillary waves in AS at least in the 
framework of a phenomenological approach. 

Now we find out the impact of surface forces on the surface entropy which 
is specific to all metalls. These are the electrostatic force and anisotropic 
pressure of electron gas in the surface zone. The contribution of electric field 
to the surface free energy is F s oc —eoE 2 5z, as discussed above. We suggest 
that the structure of expression for the free energy Eq(g) of a metal surface 
given at OK is valid for T^ but the two terms in Eq(^) are supposed to be 
temperature dependent. Then the surface excess entropy should be affected 
by the electrostatic energy and the gradient of the electron density; both 
are specific solely to the surface zone and are zero in the bulk metal. The 
standard expression S = —dF/dT yields for the entropy 

AS es = -^J^dzixw^W-eoE*) = 

(4) 
= -Sz-§ T [Xf{ d £?] + 2e,\E\5z™ 

Here we make an essential assumption that the characteristic length scale 5z 
of the surface zone remains unchanged with T whereas the squared density 
gradient and electric field should vary. This is so because the spatial density 
of surface electrons varies rapidly on the scale ~ kp as the Fermi wave- 
length ki? is T-independent. By contrast, the amplitude of the first peak 
in the electronic density p{z) is temperature sensitive, at least over T-range 



from 30 to 200 °C |lq| . According to Samojlovich's analysis [12] we associate 



the appearance of this maximum with the strength of surface electric field 



that should decrease on heating as the layering amplitude does ]15|, |ll[ . Put 
in other way, the density of over-compressed surface atoms decreases due 
to a decrease in the density gradient and electrostatic pressure, whilst the 
thickness of double layer is nearly constant. The sign of AS es depends on 
the interplay between the two terms in Eq(4) whose numerical amplitudes 
remain essentially unknown to the present study. Insofar as dE/dT is neg- 
ative AS es , may be negative too. Due to the opposite effects produced by 
the density gradient and the electric field the amplitude of AS es is expected 



to be small, hence d'y/dT is mainly affected by capillary waves especially at 
elevated T. However the theory outlined points out a possible negative sign 
in the surface excess entropy in liquid metals when AS es overcomes capillary 
wave entropy AS cap . 

3 Experimental methods 



Our light scattering technique is described in detail elsewhere J3|, |14|, 1C]. 
In brief, a beam from a 7 mW He-Ne laser (TEMoo, A=632 nm) fell on 
the liquid surface. Small-angle scattered light was optically mixed (on a 
photodetector) with a portion of the original beam, providing all the neces- 
sary conditions for optical heterodyning. The output of a photomultiplier 
was modulated at the propagation frequency of a capillary wave with the se- 
lected wavenumber q. The spectral representation of the signal was recorded 
in the frequency domain with a spectrum analyzer. The whole apparatus 
was placed on an optical table, vibration isolation being provided by four 
pressurized air cylinders in the legs. 

The liquid Ga surface was prepared in a high vacuum chamber (10~ 7 to 
10~ 8 Torr) by dropping Ga from a glass cylinder terminating in a capillary 
of diameter less than 1 mm. Liquid Ga (purity 99.999, Alfa) was dropped 
in vacuum to a Ti trough having a diameter 50 mm. The low melting rate 
ensures that a thick film of oxide remains on glass walls while pure metal 
flows through the capillary. Usually it takes some hours to complete the 
layer with strong curvature near the walls and a flat area in the middle 
part. Though no visible traces of contamination on the liquid surface were 
detected by eye the laser beam displayed a halo caused by scattering from 
a transparent oxide film which had flowed through the capillary despite all 
precautions. This continuous oxide film was fired with an Ar-ion sputter 
gun having the beam energy of 2keV at 10 -5 Torr pressure of Argon. Al- 
though the ion beam hits the area smaller than the whole surface, the film 
particles flow along the gradient of surface tension to the beam footprint. 
Hence a self cleaning mechanism due to Marangoni flow works on the liquid 
surface sputtered in vacuum. The footprint of a laser beam on the sput- 
tered gallium surface was hardly visible indicating the mirror finish and the 
absence of particles on the liquid surface (apart from an island of oxide par- 
ticles sticking to the wall). Note that the experimental w(T) dependencies 
reported below were obtained with the sputter gun on. Switching off the 
ion source caused an increase in the visibity of the laser beam footprint due 
to enhanced scattering from an oxide layer formed either due to oxidation 
by the residual gas or due to re-distribution of oxide particles on the Ga 
surface. 

Due to very high surface tension and non-wetting between gallium and 
the container a Ga-layer, usually 6-8 mm thick in the middle, easily supports 



the propagation of vibrations in the liquid sample. Although the whole 
chamber was placed on an optical table, mechanical disturbances from a 
turbo-molecular pump caused such high vibrations of the liquid sample that 
no wave spectra could be measured. We tried to obtain wetting between 
liquid Ga and the container in order to reduce the liquid depth. Wetting 
can be easily obtained by heating the sample in vacuum above 500 C. At 
this temperature Ga spreads over the container surface and thin layers can 
be prepared. However the container material very likely dissolves in gallium 
at high temperature as follows from our experiments. To keep the liquid 
sample clean one should try to achieve wetting at room temperature. 

In general, liquid gallium dropped in vacuum wets a polished container 
wall in some places but it spreads neither along the wall nor along the 
bottom. Therefore a Ga drop should be rolled along the walls until a ring of 
liquid is formed. Then the layer could be completed by shaking the liquid 
sample. After the layer is completed its form is determined by the attraction 
potential between Ga and walls; usually the surface profile was a highly non- 
symmetric (along azimuth) with concave and convex parts. Only partial 
wetting could be achieved: places that were wetted had a strong meniscus 
along the wall upwards; where no wetting took place the liquid edge was 
directed downwards. The Ga-layer prepared in such a way is about 3 to 4 
mm thick in the middle. Due to a reduction in amplitude of vibrations the 
wave spectra could be collected on such samples; wetting helps to achive 
shallow and tense layers. 

Capillary waves scatter light mainly at small angles about the reflected 
beam. The spectrum of the scattered light is the power spectrum of capillary 
waves, which is approximately Lorentzian H. The data were fitted with a 
theoretical function that incorporates the effects of instrumental broadening 
Jl7| , |i"8| ] . The spread 5q in the wavenumbers gives a corresponding broadening 
Auj in the spectrum. For the Gaussian beam the instrumental function is 
also a Gaussian 17, 18 1. The convolution of an ideal Lorentzian and the 



Gaussian instrumental function of width j3 yields [18 1 



V ; 7-oo T 2 + (a/ - co ) 2 

This integral can be evaluated in terms of the complementary error function 
of the complex argument [O, [HJ : 

S(w) = A3t[exp(-(iT/p + (u> - ^o)//5) 2 ) 

(6) 
erfc(-i[ir//3 + (a;-a;o)//3])]+-B, 

where A is the scaling amplitude and B the background. Thus, five prop- 
erties were extracted from the fit of experimental spectra: frequency loq, 
damping constant T, instrumental width /3, amplitude A, and background 
B. In the present context we concentrate only on peak frequencies uq. 
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4 Results 

The temperature behavior of peak frequencies ujq of capillary waves at the 
free gallium surface is shown in Fig. 1. To first order, the roots of the 
dispersion relation describing the propagation of a capillary wave with a 
particular wavenumber q are || 

^o = \liq z /p, (7) 

r = 2 V q 2 /p, (8) 

where 7, 77, and p are the surface tension, bulk viscosity and density, respec- 
tively. Equation (0) serves as a basis for evaluation of the tension. The data 
in Fig. 1 are fit by assuming a linear temperature dependence of the surface 
tension, 7 = a(T — To) + 70, where To is the melting point: 

w (T) = y/(a(T - T ) + 7 o)g 3 /P- (9) 

The best-fit estimate of the slope of the temperature dependence in the 
temperature interval from 100 °C to 160 °C is: a = d-y/dT = -0.57±0.09 
(mN/m)/K. The best estimate of 70 corresponding to the tension at the 
melting point (30° C) is 7o=766 ±17 mN/m. 

In the temperature range from 30 to 100 °C the temperature depen- 
dence is anomalous and the best estimate of the slope is a = 4.18 ± 0.2 
(mN/m)/K and the tension at the melting point is 7o=690±6 mN/m. The 
two data sets at temperatures above 100 C and below it were collected on the 
same sample but at two different q (295 and 310 cm -1 ) which was dictated 
by the necessity of optical adjustment of the reference beam. Theoretically 
there is no limitation for a negative slope (at T far away from the critical 
point) but there is one for a positive a. The point is that a negative ex- 
cess entropy indicates that the entropy density at the surface is less than 
that density in the bulk. We do not know how much this density changes; 
however we can estimate a maximum variation of the entropy whose spatial 
profile cannot fall below zero in the surface zone (for graphic interpretation 
see JI]). For liquid Ga at 303 K we find the bulk entropy per mole S= 
9.82+4.4=14.2 cal/mol/K §(J. Then the entropy density is S b = Sp/M Ga = 
59.7x6.1/69.7=5 MJ/m 3 /K, where MQ a is the molar weight of Ga. In or- 
der to obtain the surface entropy density this value should by multiplied by 
the thickness of surface zone: (S s ) max = S\,5z = 5 x 10 6 x 10 x 10~ 10 =5 
mJ/m 2 /K where Sz ~10 A was suggested to be the characteristic length 
scale of oscillations of the atomic density Q . This value confirms that the 
variations of the entropy occur in a relatively thick surface zone of about 
10 A indicating that the electrostatic interactions penetrate deeper than 
the topmost atomic layer. a= 4.18 mN/m/K extracted from experimental 
data at q=295 cm -1 is in accord with a =3.87 mN/m/K obtained from 



peak frequencies at 456 cm -1 on the surface of another sample (Fig. 1 (b)) 
demonstrating an increase in surface tension on heating as a real fact and 
not due to vibrational changes in uj. The damping constants inferred from 
the fit were strongly influenced by vibrations: the disperisons of V were up 
to 100 percent of the mean values. Therefore no quantitative predictions 
of the magnitude and the temperature behavior of T could be made in the 
present experiment. 

The j(T) behavior can discussed now. The effect of capillary waves on 
AS is positive and A5 cap is very likely T independent. On the other hand 
AS es becomes negative according to Eq(Q) when the condition eqE | fy |> 
Xw -§f[\{-^) 2 ] is fulfilled. This means that both the density gradient and 
the electric field decrease on heating (because a part of the surface electrons 
is injected into the bulk) setting up d^jdT < untill [ AS e i [> AS cap is 
valid. Of course, this mechanism drives the magnitude of AS untill the 
Fermi distribution of electrons is T independent; for metals it is valid up to 
1000 K. 

The light scattering experiment does not directly yield the surface ten- 
sion; only the peak frequency is directly measurable. In order to obtain 
additional information on the sign of dj/dT, the following experiment was 
carried out. A temperature gradient imposed along the sample causes a 
liquid metal flow (usually called thermocapillary convection) driven by the 
gradient of the surface tension. The surface tension gradient is indirectly re- 
lated to the liquid velocity v x via the boundary condition for the tangential 
stress at the free surface [plj[ : 

dj dv x 

dT VT = ""& (10) 

The direction of flow on the liquid surface is expected to be from hot to 
cold if the temperature derivative of the tension is negative. In the opposite 
case, the surface flow should be directed toward the hot end. Gallium was 
solidified and transferred to another vacuum enclosure having better con- 
ditions for creating a temperature gradient along the water cooled flange. 
This chamber was not an UHV: a vacuum of only about 5xl0 -6 Torr could 
be achieved. 

After sputtering the liquid surface we observed the surface flow always 
to be directed from cold to hot end in the middle of the surface. The flow 
was visible by means of light scattered from the beam footprint on the sur- 
face. Due to the extremely small scattering cross-section such observations 
are only possible when oxide particles occasionally flow through the beam 
footprint. This explaines that 7 increases with T in the temperature range 
15 to 35 °C. Further heating caused deterioration of the vacuum with the 
subsequent oxidation of the sample and quick suppression of convective flow. 
Thermocapillary convection usually couples to buoyancy convection; there- 
fore this reverse direction of the surface flow was observed only on shallow 



liquid samples. For deeper layers buoyancy convection takes over and sur- 
face flow is directed from hot to cold regardless of the temperature behavior 
of the tension. 

In sum, we studied the free Ga surface using light scattering from cap- 
illary waves. Peak frequency of waves exhibits the temperature-dependent 
behavior u(T) that corresponds to an increase in surface tension up to 100 
°C and a decrease in 7 upon heating to 160 °C. This phenomenon clearly 
indicates that the surface excess entropy switches from negative to positive 
at temperature near 100 °C. The present study shows a clear difference in 
the temperature behavior of tension between liquid gallium and mercury. 
d'j/dT of mercury, studied with light scattering technique, is continuously 
negative in the T-range from to 100 °C that is probably due to a higher 
effect of capillary waves on the surface excess entropy in mercury compared 
to the same contribution in liquid gallium. 
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Captions 

Fig. 1 Temperature variations of peak frequencies of capillary waves of dif- 
ferent wavenumbers given at each plot in cm" 1 at the free surface of 
gallium. The errors in ujq are less than the size of the data points. 
The solid lines are the best-fit solutions in the form of Eq. (||). The 
experimental points in Fig 1 (b) were collected on heating (+) and 
cooling (x). 
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